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Assuming the existence of a scalar field which undergoes "ghost condensation" and which has a 
suitably chosen potential, it is possible to obtain a non-singular bouncing cosmology in the presence 
of regular matter and radiation. The potential for the ghost condensate field can be chosen such 
that the cosmological bounce is stable against the presence of anisotropic stress. Cosmological fluc- 
tuations on long wavelengths relevant to current cosmological observations pass through the bounce 
unaffected by the new physics which yields the bounce. Thus, this model allows for the realiza- 
tion of the "matter bounce" scenario, an alternative to inflationary cosmology for the generation 
of the observed primordial fluctuations in which the inhomogeneities originate as quantum vacuum 
perturbations which exit the Hubble radius in the matter-dominated phase of contraction. 

PACS numbers: 

o 

(N 

(D 

t/J ' The inflationary universe scenario [l|, |2| is the current paradigm for early universe cosmology. While phenomeno- 

logically successful, inflationary models face various conceptual problems (see e.g. Q), most notably the singularity 
problem J4|. Hence, it has been of great interest to construct nonsingular bouncing cosmologies (for a recent review 
see e.g. Q). 

It is possible that a nonsingular cosmological bounce is followed by a period of inflation (see e.g. p for an example). 

f-^ , However, it is also of interest to explore alternatives to inflation which could result from a bouncing cosmology. In 

fli ' 0) @] it was realized that quantum vacuum fluctuations which exit the Hubble radius in a matter-dominated phase of 

r^ , contraction obtain a scale-invariant spectrum of curvature fluctuations on super-Hubble scales. Provided that on the 

'~~' • large scales corresponding to current cosmological observations the fluctuations pass through the bounce unaffected by 

-v^ \ the new physics which yields the cosmological bounce, then we have a mechanism alternative to inflation for producing 

►^ . a scale-invariant spectrum of cosmological perturbations. This alternative is called the "matter bounce scenario" (see 

^^ ' e.g. [9] for introductory reviews). The matter bounce scenario can be distinguished from simple inflationary models 

lO • through the shape and amplitude of the bispectrum, the three-point function of cosmological perturbations, which 

^^D ' has a special shape and an amplitude which is large compared to what results in simple inflationary models [10|. 

^N \ From the Penrose-Hawking singularity theorems [11[ it is clear that new physics is required in order to obtain 

t"~~. ' a cosmological bounce. There have to be modifications either to the gravitational action, or else matter which 

^D ' violates the usual energy conditions has to be introduced. There have been various suggestions on how to obtain a 

^^ I cosmological bounce by means of modifying gravity. For example, the "non-singular Universe" construction of fl2| 

", , ' yields a cosmological bounce, as does the ghost- free higher derivative gravity action of [1^ . It is also possible to obtain 

^ , a bounce in the context of "mirage cosmology" where our space-time is a brane moving in a curved higher-dimensional 

background space [ij . As was realized in [l^ , in the presence of spatial curvature a cosmological bounce also occurs 

^ . , generically in Hofava-Lifshitz gravity [16|. 

j^ ■ If we maintain General Relativity as the theory of space-time, then a cosmological bounce can also occur if we 

suitably modify the matter sector. A generic way of obtaining a bounce is to add to the usual matter sector which obeys 

the standard energy conditions a "ghost" sector with opposite sign kinetic ener gy t erms in the matter Lagrangian. 

This class of models is called "quintom" models \\T\ and can lead to a bounce [l^. One way to realize a quintom 

bounce is by considering a scalar field ^ with standard kinetic term and with mass m and adding a ghost scalar field 

(^ with mass M ^ m. This model was considered in detail in |19l|. A specific realization is in the context of the scalar 

field sector of the Lee- Wick Standard Model !20l, as studied in [2l| |. A serious challenge for all models with a ghost 

field is the ghost instability problem (see e.g. 22]). 

All analyses of the evolution of cosmological fiuctuations thr oug h non-singular bounces (e.g. [2j| in the case of the 
model of [I^l , [ij| in the case of a mirage cosmology bounce, [I^ [21| in the case of quintom models and [2J] in the 
case of the Hofava-Lifshitz bounce) indicate that the spectrum of the curvature fluctuations does not change during 
the bounce phase on length scales larger than the duration of the bounce. 



There are two serious problems for matter bounce scenarios. The first is the fact that some bounces are unstable 
to the presence of radiation, the reason being that the energy density of radiation increases faster as the universe 
contracts compared to the increase in the effective energy in the terms yielding the bounce. In particular, this was 
shown [25J to be a serious problem in the case of the bounce in the Lee- Wick model, and the problem will likely also 
be serious in other models in which the ghost field yielding the bounce is a scalar field. In contrast, models which are 
based on asymptotically free gravity models such as [I4I and [13[ will be free from this problem, as is the bounce in 
Hofava-Lifshitz gravity. 

An even more serious problem is the instability of the contracting pre-bounce phase to the presence of anisotropics 
since the energy density in anisotropies scales as a~^ (a(0 being the cosmological scale factor) compared to a~^ for 
radiation ^ . This implies that the cosmological bounce cannot be described in terms of a homogeneous and isotropic 
background. 

In this paper, we present a ghost condensation model in which a cosmological bounce can be obtained. Since there 
is no ghost in the perturbative spectrum of the theory, the model is free from the ghost instability problem of [24I . 
The model is also free from the two problems we just mentioned above. The effective equation of state of the ghost 
condensate is w > 1 and hence, as studied in ^2T\ in the context of the Ekpyrotic scenario [285, the anisotropies do not 
come to dominate and the cosmological background remains well described in terms of a homogeneous and isotropic 
metric. 

Ghost condensation was proposed [23 as a theoretically consistent modification of gravity in the infrared. It is 
based on the analog of the Higgs mechanism in the kinetic sector of a scalar field Lagrangian. The kinetic term X of 
the scalar field cj), i.e. 

X = -5^''9^0a,0 (1.1) 

where we are choosing the metric signature to be (-, + + +), appears in the Lagrangian in terms of a function P{X) 
with a nontrivial minimum, such that the leading term in X in the Lagrangian has the "wrong" sign, i.e. appears as a 
ghost. If the function P{X) has a non-trivial minimum at X = c^ (the "ghost condensate") and the field configuration 
takes on this value, then fluctuations about this configuration have normal kinetic term. Thus, there is no ghost in 
the spectrum of fluctuations about the ghost condensate. 

It was soon realized that, although initially introduced [29| to provide a new explanation for infrared effects in 
gravity, the ghost-condensation mechanism has many applications in cosmology at ultraviolet scales. For example, 
it can be used [3(| to provide a new mechanism for inflation. It was realized in [3l| that the ghost condensate 
mechanism can also provide stable violations of the null energy condition and thus used to construct non-standard 
cosmologies, including bouncing ones. Ghost condensation was used in [32, l33J to provide non-singular versions of 
Ekpyrotic cosmology. In this paper we make use of the ghost condensation mechanism to construct a non-singular 
matter bounce. 

The outline of the paper is as follows. In the following section we briefly review the idea of ghost condensation. In 
Section 3 we work out the requirements for the ghost fleld potential V{4>) in order to obtain a bounce, in particular a 
bounce which is stable against the addition of radiation and anisotropic stress. In Section 4 we discuss the transfer of 
cosmological fluctuations through the ghost bounce phase. The final section contains our conclusions and a discussion 
of the results. 

To set the notation, we work in terms of a spatially fiat Friedmann-Robertson- Walker background with metric 

ds^ = -df + a{tf{dx^ + dy"^ + dz^) , (1.2) 

where a(t) is the scale factor, t is the physical time and x, y and z are the co-moving spatial coordinates. The Hubble 
expansion rate is H(t) = a/a. 

II. GHOST CONDENSATION 

The Lagrangian of ghost condensation takes the following general form 

C = M*P(A:) - 1/(0) , (2.1) 

where M is the characteristic mass scale of ghost condensation and we follow the usual convention in the ghost 
condensation literature to take (j) to have dimensions of length such that X is dimensionless. To obtain a ghost 



^ This leads to the BKL l26l chaotic behavior at the "Big Crunch" singularity. 



condensate, the function P{X) must have a non-trivial minimum. To avoid a cosmological constant problem, we set 
the value of P at the minimum to be zero. We will use the prototypical example 

P{X) = ^{X - c')\ (2.2) 

where c is a dimensionless constant which has nothing to do with the speed of light (we are following the notation 
of the original papers on ghost condensation). The ghost condensate necessarily breaks Lorentz invariance. The 
homogeneous and isotropic ground state corresponds to the background scalar field configuration 

(j) ^ ct. (2.3) 

The equation of motion of the ghost condensation field is 

IP' + 202 p"U + 2,HP'^ = -M""— - , (2.4) 

where P' = fy- 

The analogs of the Friedmann equations obtained by coupling the homogeneous and isotropic ghost condensate to 
Einstein gravity are 

iM^H^ ^ M'^{2XP' -P)+V + p^, (2.5) 

and 

2MlH = -2M^XP' - (1 + w^)pra . (2.6) 

We have assumed that in addition to the ghost condensate, there is also regular matter with energy density pm and 
equation of state w — w,n , where w is the ratio of pressure density to energy density. 
Let us now consider homogeneous fiuctuations tt about the ghost condensate field, i.e. 

<^(i) = ct + TT{t). (2.7) 

Inserting into (J2.5I) and (12. 6p we obtain the following expressions for the energy and pressure densities of the homo- 
geneous TT field: 



px = M^^tt{1 + 0[-))+V, (2.8) 

c 

Px+Px = MV7t(1-HO(-)). (2.9) 



In order to obtain a cosmological bounce, it is necessary to cancel the positive energy density of regular matter with 
the negative energy density of the ghost condensate. This can be achieved by having negative tt. To complete the 
system of equations for the fiuctuations of the ghost condensate, we write down the variational equation with respect 
to <^. To leading order in 0{ — ) it is 

c^a-^'dtia^i:) ^ -2M-^^. (2.10) 

a<p 

By inserting the ansatz (|2.7p into the Lagrangian (j2.ip it follows that P' + 2(p^P" ^ is the necessary condition 
for the theory expanded about the condensate to be ghost free. 

III. BOUNCE INDUCED BY A GHOST CONDENSATE 

In the absence of regular matter, a ghost condensate leads to a cosmological bounce if the ghost field potential 
energy (which we assume to be positive) is cancelled by the negative ghost condensate energy. To leading order in 
0{-) the condition for the bounce point is 

mVtt = -V . (3.1) 



The second Friedmann equation, i.e. (j2.6p . then imphes that H > a,t the bounce point, i.e. the scale factor indeed 
makes a transition from a contracting phase to an expanding phase. 

To obtain a bounce, it is necessary to have a non-trivial potential V{(l)) since in the absence of a potential the 
equation of motion (j2.10p implies 

a^TT — const (3-2) 

and hence the bounce condition, which in the absence of a potential is tt = 0, cannot be reached. 

To obtain a matter bounce, we must consider the equations of motion in the presence of regular matter. We assume 
that the universe begins at very large negative times with the ghost condensate in its ground state X = c^ (or tt = 0). 
The energy density in regular matter increases in time as 

p„(t) ^ a(i)-3(i+"'"' . (3.3) 

We are interested in the cases of cold matter (w„j — 0) and radiation {wm — 1/3), or in anisotropic stress for which 
Wm — 1- From (|3.2p it follows that in the absence of a potential V for (p 

px - a-^ . (3.4) 

Hence, the ghost energy density cannot catch up with the matter energy density and no bounce is possible. Thus we 
conclude that, as in the case without matter, also in the presence of matter it is necessary to introduce a potential 
V{(j)) in order to obtain a bounce. 

In order to obtain a bounce, the energy density in the ghost field must increase sufficiently fast: 

PX ^ a(t)-^ (3.5) 

where p > 3 is the minimal requirement, p > 4 is required if the bounce is to be stable against the presence of regular 
relativistic radiation, and p > 6 is required if the bounce is to be stable against the presence of anisotropics. 
In order for the ghost energy density to increase faster than a"'^, we require 

dt{a^Tr) < (3.6) 

and hence from (j2.10p 

dV 

^>0. (3.7) 

We assume that the ghost field (p starts at large negative values. Hence, we make the following ansatz for our 
potential 

F(0) - VoM-"cl>-~" , (3.8) 

where Vq is a constant with units of potential energy density which sets the overall scale of the ghost energy density, 
and a is a constant to be determined by the requirement that the ghost energy density increases sufficiently fast. 
This potential diverges at (/> = 0. However, we know that the ghost condensate Lagrangian should only be trusted at 
energy scales lower than M ^. Hence, we will assume that the divergence of the potential is cut off when the potential 
reaches this value. For example, we can assume that the potential tends smoothly to the limiting value M^ which is 
taken on at (^ = (see Fig. 1 for a sketch of this potential). 

Let us now derive the condition on a for which a non-singular bounce which is stable even in the presence of 
anisotropics will arise. For this, we return to the equation of motion (J2.10D which after inserting the ghost fluctuation 
ansatz (j2.7p and expanding to leading order in n becomes 

c^dtia^jr) = -2aHr^'"^ . (3.9) 

a<p 

Inserting the form of the potential (j3.8p . the ansatz for the field from (|2.7p . and working in the range of times for 
which \ct\ ^ \-k\ yields 

n + 'dHir = 2c-2vbM-*-"a(ci)-("+i) . (3.10) 




FIG. 1: Sketch of the effective potential of the ghost condensate field. 

We want to derive the condition on a such that in a background dominated by matter with an equation of state 
parameter w, the energy density of the ghost kinetic term which is proportional to tt increases faster than the energy 
density in matter as the universe contracts. Hence, we insert the following form of the scale factor 

a(t) = (^)2/(3{i+-)) (3.11) 

into p.lOp . We see that the source term (the right hand side of (I3.10p leads to 

TT - t~". (3.12) 

and the condition for the ghost kinetic energy density to grow faster than the matter energy density is 

a > 2. (3.13) 

Since we wish the potential to be an even function of 0, we can choose a — A. In this case the model is marginally 
stable against the addition of anisotropic stress. The case a = 6 would provide a model which is stable. 

Let us now summarize the evolution of the background cosmology in our ghost bounce model. We begin at t — >■ — oo 
in a contracting, matter-dominated phase with ghost condensate field — >■ — cxd and the ghost condensate in its ground 
state with TT = 0. The energy density of the universe is initially dominated by non-relativistic matter. As the universe 
contracts, first radiation with w = 1/3 and then anisotropic stress with w = 1 take over. However, the energy density 
in the ghost condensate grows faster than any of these energy densities and, at some point sufficiently close to t — 
begins to dominate. Note that both the negative energy of the ghost kinetic term and the positive potential energy 
of the ghost field increase at the same rate. Eventually the ghost field (j) reaches the value when the potential fiattens 
out. At this point, the negative kinetic contribution to the ghost energy catches up, and the bounce point H = 
is reached. Since at this point tt < we have H > and the universe begins to expand. The ghost field evolves 
to positive values, and both the ghost potential and kinetic energies begin to decrease rapidly, allowing the energy 
densities of anisotropic stress, radiation and cold matter to dominate again. 

The evolution of the background has been studied numerically and the results are presented in Figures 2-5. In 
all four figures, the horizontal axis is time in Planck units. The vertical axis is 4>, (j), and H respectively. There is a 
nonsingular bounce at time i = 0. The graphs are based on solving Eqs. (12. 4p and (|2.5I) using Mathematica, beginning 
the evolution at the bounce point i = where "initial" conditions (again in Planck units) 0(0) = and 0(0) = ■\/2/3 
are chosen. The value of M was taken to be M = 2 x 10"'^, and the potential was taken to correspond to a = 4, with 
a value of Vq such that 

n0) = ^r^. (3.14) 

The energy density of matter sector was taken to be pm(0) — 10^^^. As is obvious from the graphs, a smooth 
nonsingular bounce results. 
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FIG. 2: Evolution of the ghost field (j) as a, function of time in our numerical simulation. 



-20000 -10000 




10000 20000 



FIG. 3: Evolution of the ghost field velocity as a function of time in our numerical simulation. 



IV. A FIRST LOOK AT COSMOLOGICAL PERTURBATIONS IN THE GHOST CONDENSATE 

MODEL 

A. Generating a Scale-Invariant Spectrum 

In this section, we will show that the perturbations induced by the ghost condensate do not change the primordial 
spectrum of curvature fluctuations induced by regular matter in the contracting phase. Thus, in the matter bounce 
background the scale-invariant spectrum of fluctuations on scales which exit the Hubble radius during the mater- 
dominated phase of contraction will not be distorted by the ghost condensate, neither in the contracting nor in the 
bounce phase. First, in this subsection, we review how a scale-invariant spectrum of curvature fluctuations arises in 
the matter bounce scenario. 

To set the framework, we present in Figure 6 a space-time sketch of a matter bounce. The horizontal axis is space, 
the vertical axis is time, and the bounce time is t = 0. We are interested in fluctuations on scales which are currently 
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FIG. 4: Evolution of the Hubble parameter H as a. function of time in our numerical simulation. 



observable, and assume that these scales exit the Hubble radius in a matter-dominated phase of contraction. This 
assumption is not very restrictive. If the bounce cosmology background is time-symmetric and the expanding phase 
corresponds to our currently observed universe, then the assumption is satisfied for all scales on which the primordial 
matter power spectrum is currently well measured ^. If the bounce is not time-symmetric, then according to the 
Second Law of Thermodynamics the expanding phase will have more entropy than the initial collapsing one and the 
radiation phase is shorter than in the case of the symmetric bounce, which means that a wider range of scales exits 
the Hubble radius in the matter-dominated phase. 

As shown in the Figure, scales of cosmological interest today exit the Hubble radius during the matter-dominated 
phase of contraction. If the bounce takes place at an energy scale comparable to that of Grand Unification, then the 
current Hubble radius corresponds to a physical length of about 1mm at the bounce. This scale is in the far infrared 
compared to the scale which sets the physics of the bounce. 
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FIG. 5: A sketch of the evolution of perturbations with different conioving wave numbers k in the matter bounce. The horizontal 
axis represents comoving coordinates, the vertical axis is time, with the bounce point being t = ts- The two vertical lines 
correspond to the wavelengths of two different scales, the thick black curve gives the Hubble radius. 

The theory of cosmological perturbations in the model with both regular matter and a ghost condensate in general 
FRW background was worked out in [Sa] ■ Our goal here is to show that the presence of the ghost condensate and its 
induced curvature fluctuations will not mess up the scale-invariant spectrum of perturbations resulting from regular 
matter in a matter-dominated phase of contraction. We have to show first that the ghost fluctuations in the matter- 
dominated phase of contraction do not grow faster than the fluctuations without the presence of the ghost field. 
Secondly, we need to show that the curvature fiuctuations which are induced by the coupling between regular matter 
and the ghost field are negligible, and thirdly, we need to show that the spectrum of fiuctuations on scales of relevance 
to current cosmological observations is not changed during the short phase where the ghost condensate dominated 
the dynamics. 

We begin with the general framework of metric fluctuations. Working in in longitudinal gauge (see [36] for a 
comprehensive review of the theory of cosmological perturbations) the metric including scalar fluctuations take the 
form 



ds' 



(1 + 2$)rff2 + a(f)2(l - 2*)dx2 



(4.1) 



where <I> and ^ are the two gravitational potentials which depend on space and time and encode the metric pertur- 
bations. Since there is no anisotropic stress at linear order (both in the regular matter sector and in the ghost matter 
Lagrangian), we have ^ = $ as a consequence of the perturbed off-diagonal spatial Einstein equations. The remaining 
potential $ is the relativistic generalization of the Newtonian gravitational potential. 

To describe fluctuations on super-Hubble scales, it is more convenient to use the fluctuation variable v in terms of 
which the action for canonical fluctuations has canonical form [371]. If matter is a simple scalar fleld ip, then v can be 



^ This background predicts a kink in tlic matter power spectrum on scales which can be probed by Lyman a observations - see [SJ 
discussion of this point. 



for a 



expressed in terms of the matter and metric fluctuations 6lp and $ via 

V = a[5Lp+-^] (4.2) 

where z is a function of matter and metric background whose form does not concern us here (see [36[ or [38| for an 
introductory review). Closely associated with v is the variable ^, the curvature fluctuation in co- moving gauge (called 
unitary gauge in the ghost condensate literature): 

C = a-^v. (4.3) 

In an expanding universe the dominant of the two modes of C has constant amplitude. However, in a contracting 
universe the dominant mode of C, is increasing on super-Hubble scales. In fact, in a matter-dominated phase of 
contraction the dominant mode of v scales as (see [3, Q for the original works and Q for recent reviews) 

v{t) ^ i-i/3 , (4.4) 

and hence 

m - 1-' ■ (4.5) 

The above is exactly the growth rate needed to convert an initial vacuum spectrum for fluctuations on sub-Hubble 
scales into a scale-invariant spectrum on super-Hubble scales. We briefly review the argument. The vacuum spectrum 
for the canonical fluctuation variable states that 

v{k,t) - fc-i/2 (4.6) 

on sub-Hubble scales in the contracting phase. The condition for the time tnik) when the scale k crosses the Hubble 
radius is 

tH{k) = k-^a{tH{k)). (4.7) 



Making use of the scaling (j4.5p , the vacuum initial conditions (j4.6p which hold until Hubble radius crossing, and the 
Hubble radius crossing condition (|4.7p we obtain 

C(fc,i) = ^-^a{tH{k))v{k,tHik)) ^ tHiky/'vitnik)) ^ k-^k-''\ (4.8) 

which corresponds to a scale- invariant power spectrum. Note that in the last step we have used the scaling tnik) ^ k~^ 
which holds in a matter-dominated phase of contraction. 

Thus, we have shown that vacuum fluctuations which exit the Hubble radius during a matter-dominated phase of 
contraction acquire a scale-invariant spectrum of perturbations on super-Hubble scales. This spectrum is maintained 
if the background equation of state changes when radiation starts to dominate. We will have to show that the presence 
of a ghost condensate does not change the evolution of the fluctuations in the contracting phase, and that the spectrum 
is not distorted in the phase dominated by a ghost condensate. 

B. The Evolution of Ghost Perturbations 

Let us now turn to the ghost matter sector. Introducing a fluctuating fleld 7r(i, x) which depends both on space 
and time, the Lagrangian to quadratic order in tt about a non-trivial background ttq (the -k appearing in the previous 
section) is [3l| 



lM^ir-+HMl,{V.f-\r 



C = -M^^ + HM^iiVirf - -M^iV^irf , (4.9) 



where V is the spatial gradient operator. Note that to leading order in tt, no quartic order spatial gradient terms 
app ear if we only use the Lagrangian (12. ip . and to obtain (14.91) we had to introduce higher derivative terms, following 
[23,l31|. To derive the form of the second term in (|4.9p . we made use of the ghost background equations of the previous 
section. 



By working in unitary gauge in which ( appears directly as the fluctuation of the diagonal term in the spatial 
metric, it was shown in Section 3.3 of |31| that for long wavelength modes 

C r^ a-^ r^ t-^ , (4.10) 

where we have used the matter-dominated scaling of a{t) in the last step. By comparing (j4.5p and (I4.10p we see that 
the matter and ghost-induced fluctuations grow at the same rate on large scales. 

To show that the presence of a ghost condensate sector will not change the dominant contribution to the spectrum of 
cosmological fluctuations in the phase before the ghost condensate begins to dominate, we need to study the spectrum 
of the metric fluctuations induced by the presence of coupling between the matter and ghost sector. To do this, we 
make use of the equations for the metric fluctuation variable $ which were derived in [35| . 

As matter we consider a ghost condensate in addition to a matter fluid, and we focus on a matter-dominated 
contracting background. Thus, the Einstein equations are 

M^,G,.^TlfJ+T,,, (4.11) 

where we use a superscript {(f>) to denote the ghost condensate part. We split $ into the contribution $,„ which would 
be obtained in the absence of matter fluid, and a term $g which is induced by the presence of the ghost: 

* = *™ + $<, . (4.12) 

If we neglect the anisotropic stress and entropy perturbation, the general equations of motion for $ are [35| 



dt^g+H^g = 



a k 

.M2 2 k2 



d.X= ij^-J^^)'^. + S,^ (4.13) 



where a = M"^ /M"^ and the source term is given by 



'pi 






-Z2^r. = ^ . (4.14) 



Note that the two equations in (|4.13p come from combining all of the perturbed Einstein equations, i.e. the — 0, 
— z, and i ~ i components. Thus, the set of these equations is complete. Combining these two equations we obtain 
the following equation of motion for the Newtonian potential: 

This equation can be simplified by factoring out the expansion/contraction of space via defining the new variable 
$ = a{t)^, and by working with conformal time r defined via dt — a{t)dT. Then, the equation of motion for the 
Newtonian potential can be rewritten as 

where a prime indicates the derivative with respect to t. We have already used the matter contracting background 

a — J- (where to is some reference time) to get the above equation. 

We will solve (I4.16P in the large wavelength approximation in which we neglect the k'^ term. Then, the general 
solution is a combination of the general solution ^^ of the homogeneous equation and a particular solution $^ of 
the inhomogeneous equation which can be determined by the Green function method. To simplify the notation we 
introduce the constants where 

^ = ^ 
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The general solution of the homogeneous equation is 

1.^ = cie^Vir j^ cae-'^v^^ , (4.18) 

where ci and C2 are constants which multiply the two basis solutions of the equation Ci(t) and C2(t), respectively. 
The Wronskian e(r) obtained from the two basis solutions is 

e(r) = {2k^r^ (4.19) 

in terms of which the particular solution $p becomes 

$P(t) = Cl{T)JdT'e{T')C2{T')S{T')~C2{T)fdT'e{T')Cl{T')S{T'), (4.20) 

where S{t) is the source term in (J4.16I) . Recall that the general solution for $„ on super Hubble scales j8i] is 

$,„ = a(<)$™ = Dt^ + ^ , (4.21) 

where D and S are constants. Inserting this equation into the source term in (j4.16p . and approximating the mode 
functions Ci(r) and C2(t) of the homogeneous equations as constant (which is justified over short time intervals for 
long wavelength modes), we obtain the following approximate form for the particular solution 

$^(r) c. fc^7r(f r^ - f r-^) . (4.22) 

Note that the effects of the fc^ term in the equation of motion (|4.16l) can be included to leading order using the Born 
approximation. In this approximation, the contribution to "I>g induced by the fc** term is obtained by taking the fc'^ 
term to the right-hand side of (|4.16p . evaluating it for the homogeneous solution $g, and using it as a second source 
term in (|4.20p to obtain a second contribution to the particular solution. This source term scales as r""*, and hence its 
contribution to <I>^ will scale as t~^, it also grows slower than $„. On the other hand, the spectrum is suppressed by 
a factor of k^ and is thus highly blue and completely negligible on scales relevant to current cosmological observations. 

Returning to the contribution to $p which dominates on large scales, we conclude that on super-Hubble scales the 

leading term comes from ^p and is 



Comparing with the matter perturbation 



1,2 c 

$, ^ a(i)-i$, ^ ^. (4.23) 



*m ^ 4 ' (4-24) 



we see that the perturbation spectrum of the terms generated by the ghost condensate is blue, and that it grows more 
slowly than the contribution from the matter sector. This is a self consistent result. As analyzed at the beginning 
of this section, the spectrum of the perturbations induced by matter is scale-invariant because the amplitude of 
the perturbations which cross the Hubble radius at an earlier time grows and catches up to the amplitude of the 
perturbations which cross the Hubble radius at a later time in the contracting phase. Since the perturbation induced 
by the ghost condensate grows slower than that induced by matter, the perturbations which cross the Hubble radius 
at an earlier time cannot catch up to the amplitude of the perturbations which cross the Hubble radius at a later 
time, and thus the spectrum is blue. The scale-invariant spectrum of the perturbations induced by matter will be 
preserved in the contracting phase. 

The third step is to show that the spectrum remains unchanged around the bounce phase. We only need to show 
that the amplitude of the curvature fluctuation induced by the ghost condensate will not grow too much during the 
bounce phae. This can be seen easily from Ea. (l4.15p . Since the duration of bounce phase is very short and scale 
factor a{t) is almost a constant we can parameterize the background as follows: 

H = e-{t-tB) , (4.25) 
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where is is the time at the bounce point, and is a constant with :§> H^, where He is the value of H at the matching 
surface between the contracting phase and the bounce phase ^ . On scales of interest to current observations we have 
^ <^ Hi <^9, and thus Ea. (l4.15p can be rewritten as 

92$^ + ^$^ =0, (4.26) 

the solution of which is 

$g = die*^* + d2e-'^* . (4.27) 



Since the bounce phase is so short we thus have 



$,(i) ^ $=, (4.28) 



also at the end of the bounce phase (where $g is the Newtonian potential induced by the ghost condensate at the 
beginning of the bounce phase). The perturbation induced by the ghost is almost the same as at the end of the 
contracting phase. 

Note that near the bounce point the perturbation re-enters the Hubble radius for a very brief time interval, but 
/c^ ^ 6* is still true since scale factor is almost a constant during the bounce phase, and thus Eq. (|4.26p is still valid 
for describing the perturbations near bounce point. We thus see that, as long as the perturbations induced by the 
ghost are unimportant during the contracting phase, they will also be unimportant during the bounce phase. 

Previous work in the context of a wide range of non-singular bouncing models (see e.g. [Ij, [T^ |2l|, |2j, UM) ^^^ 
shown that the matter-induced fluctuations do not change during the bounce phase. Thus, we conclude that the full 
spectrum of cosmological perturbations in our model does not change during the bounce phase. 

C. Stability of Ghost Perturbation 



As derived in Section 3.3 of 31[, the action for the curvature fluctuation variable C, is 



S ^ f d^xdta^t) U(i)C' + S(0( — )' + C{t){^f] , (4.29) 

J a a^ 

where the coefficient functions A{t),B{t) and C{t) are given in the Appendix. Near the bounce point we can take 
the scale factor to be a constant and neglect the terms proportional to H which appear in the coefficient functions. 
Then, for each Fourier mode we obtain a harmonic oscillator equation, and the dispersion relation is 

2 -JM'M* + 4M4g)P + 2M^iM^k^ 
" = 2M^^ ■ ^^-^"^ 

As is obvious from (I4.30p . there is an instability for all long wavelength modes. In fact, the instability is a 
combination of two instabilities which can be seen [29l l31| at a heuristic level, the first being a gradient instability 
which appears if 7J > as it is near the bounce, and the second a Jeans type instability. These are due to the second 
and first terms, respectively, in the coefficient of the fc^ term in (|4.30|) . It is easy to determine the value of k for which 
the instability is the strongest. The exponent lOc of the instability is 

I MAP ■ Mil 

uj^^ -——— + H ^. 4.31 

4 Mp2 M^M ^ ' 

We will now argue that the bounce period is too short for this instability to change the spectrum of cosmological 
perturbations. To estimate the time period At of the bounce, we use 

At(/i = A0, (4.32) 



Note that this implies that the duration 5t of the bounce phase satisfies \/0(5t <C 1. 
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where A(/) is the interval of corresponding to the bounce. We will estimate this distance by the field value for which 
the potential p.Sp approaches its maximal value which we take to be ¥{(/)) = Af*. This gives 

Ac^r.M-'{^S^\ (4.33) 

Inserting (p c^i c and setting the constant c = 1 we obtain 

The first term in the square bracket is clearly much smaller than one. To estimate the value of the second term, we 
make use of (j2.6p and (j2.9p to estimate the value of H: 

-^ 

with which the second term in the square bracket of (|4.34[) is of the order of tt which cannot be larger than order 
unity. Thus, if Vb < iW then 

Aiwc < 1 , (4.36) 

and the instability in the ghost condensate phase does not have time to develop. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have presented a realization of the "matter bounce" by means of a ghost condensate. Compared to 
other realizations of the matter bounce using a modification of the matter sector, our model has several advantages: 
first of all, there is no ghost in the perturbative spectrum of states. Secondly, the background cosmology is stable 
against the addition of both regular radiation and of anisotropic stress, the latter implying that our model will likely 
be free of the chaotic mixmaster behavior which plagues many models. 

We have also studied the evolution of cosmological perturbations through the bounce. We considered metric 
fluctuations on super-Hubble scales in a model with both background matter (a perfect fluid with equation of state 
p = 0) and a ghost condensate. We have shown that the ghost-induced fluctuations grow less fast during the contracting 
phase than the curvature perturbations in a pure matter model, and that in addition they have a blue spectrum. 
Hence the scale-invariance of curvature fluctuations on super-Hubble scales in the contracting phase (assuming that 
the inhomogeneities originate as quantum vacuum fluctuations) is maintained. Since the curvature perturbations 
which are induced by the gravitational coupling between ghost condensate and regular matter have a blue spectrum 
they are irrelevant on scales of interest to current cosmological observations. Finally, we have studied the evolution of 
the fluctuations through the bounce phase, a phase dominated by the ghost condensate. There is a gradient instability, 
but the bounce phase is sufficiently short such that this instability has no time to develop. 

In conclusion, we have shown that our ghost condensate background cosmology provides a realization of the matter 
bounce alternative to inflation for generating a scale-invariant spectrum of adiabatic fluctuations. 

Ackno'wledgments 

This work is supported in part by a NSERC Discovery Grant, by a grant from the FQXi, and by a Killam Research 
Fellowship awarded to R.B. We would like to thank Shinji Mukohyama, Yi Wang, and Wei Xue for useful discussions. 

Appendix A: Jeans instability during late time evolution 

In this subsection we discuss the Jeans instability during late time evolution. The scalar excitation of the ghost 
condensate gives rise to an instability of the vacuum analogous to the Jeans instability of pressureless matter coupled 
to gravity. The length and time scales associated with the Jeans instability are 

L.-% Tj^% (Al) 
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Any fluctuation with length scale larger than the Jeans length will be unstable. Jeans collapse produces lots of lumps 
of scalar excitation, and the universe will be filled with lumps. These lumps may produce phenomenona which are 
inconsistent with current cosmological observations, like excess lensing of light, supernova time delays and so on. 
These effects increase in magnitude as M increases. If matter is dominated by the ghost condensate field as will 
be the case if the ghost condensate is used to model dark energy, an upper bound on the energy scale M of ghost 
condensation has been derived in Ref. |39l |: 

M < lOOGeV. (A2) 

However, in our model in which the ghost field potential term plays an important role, the ghost condensate only 
dominates matter around the bounce point. Thus, our model is free from the above constraint. For example, if we set 
a = 4 which corresponds to the ghost condensate being marginally stable against anisotropic stress, then the energy 
density of the ghost condensate scales as a~^, and will be diluted away rapidly at late times in the post-bounce phase. 
It is a standard result in the theory of cosmological perturbations (see e.g. [40, |41| for textbook discussions) that 
fiuctuations in a subdominant component are held back and cannot grow more than logarithmically on scales where 
the fluctuations in the dominant component are stable. For example, matter fluctuations in the radiation epoch of 
Standard Cosmology only grow logarithmically on sub-horizon scales. This implies that the ghost condensate does 
not clump in the radiation phase of Standard Cosmology, and that the ghost fluctuations track the usual matter 
fluctuations in the matter epoch. Hence, in our model there are no constraints from late time cosmology on the scale 
M. 

Appendix B: Various Coefficients 

In this subsection we give the expressions for the coefficient functions A{t),B{t) and C{t) which appear in (J4.29I) . 
They are taken from Section 3.3 of J3l| : 

2 M^ (m^ - 9M^ H^-2 A/2 ^"j 

A{t) = ^ —^ ^ (Bl) 

4Af4jiJ2 + M M4 



M, 



Bit) = i^i-^^ -24Aff,M m +M i^r - 2Mi,H)(M'^M - AM$,,H) 

(4A//4ji/2 + M'Af4)2 L 

+ AH'^iM'^Ml.Jd^ + 4:]V4iH) - ^MI^M^Hh\ (B2) 



2A'f4,Af^ 



Cit) = r!''^-i . ■ (B3) 



AM^^ff^ ^ M^ M^ 
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